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We study machines that take N identical replicas of a pure qudit state as input and output a 
set of Ma clones of a given fidelity and another set of Mb clones of another fidelity. The trade-off 
between these two fidelities is investigated, and numerous examples of optimal N —* Ma + Mb 
cloning machines are exhibited using a generic method. A generalisation to more than two sets of 
clones is also discussed. Finally, an optical implementation of some such machines is proposed. This 
paper is an extended version of [xxx.arxiv.org/abs/quant-ph/0411179]. 

in ! 
o ■ 

O ' !• INTRODUCTION 

<n : 

^"V The no-cloning theorem pj states that it is in general impossible to perfectly clone the state of a quantum system. 
This no-go result is a typical feature of quantum information, and is deeply related to quantum error correction |2j, or 
the impossibility of super-luminal signalling 2j . But the no-cloning theorem can also be turned to a valuable resource 
as demonstrated by quantum cryptography [4|. 

Although perfect cloning is forbidden by quantum mechanics, it is of fundamental importance to analyse how well 
we can approximately copy a quantum state into several quantum systems in order to understand how quantum 
_^ ■ information distributes. The simplest instance of this problem is the duplication of the pure state of a qubit, and 
' has been considered in [j| (a qubit is a two- level quantum system) . Many generalisations and variants have followed. 
Among them are the issue of producingM clones of a qubit state from N input replicas Q and its generalisation 
to qudits (d-level quantum systems) 013], cloning from non identical input states 0, or non- universal cloning [ldj . 
Cloning of continuous variable systems has also been considered |TT| . 

This paper deals with universal N — ► Ma + Mb + • • ■ asymmetric cloning of pure qudit states. That is, we are 
given N identical replicas of an unknown pure qudit state, and we want to produce Ma approximate clones all with 
a same fidelity, say F A , Mb clones, all with a same fidelity F B , .... We are interested in the trade-off between the 
qualities of the various sets of clones. Such machines have already been examined in [TJ, but only 1 — > 1 + 1 cloning 
was considered, and optimality was only proven in the case of qubits. N — > Ma + Mb cloning machines have at 
least two interesting applications. First, in the case where Ma = N and Mb — > oo, we get a machine that allows to 
study the trade-off between the acquisition of knowledge about the state of a quantum system and the disturbance 
undergone by this system. Second, some N — > Ma + Mb cloning machines have been proven to be a useful tool when 
^ | investigating the security of some quantum key distribution schemes 0, . 
O" 1 Here is a summary of our results. 
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• A generic method to get optimal N — > Ma + Mb cloning machines of pure qudit states is prescribed for any 
^ ■ values of N,M A ,M B ,d, and a 'natural' conjecture of what cloning machines this method should always produce 

is proposed . This conjecture has been supported by all the cases we have examined. These results are essentially 
the content of Section [H] 

• With the aforementioned method, optimal asymmetric cloning of qubits has been investigated in the following 
cases: 1 — * 1 + n cloning of qubits (see section Ull CI) and n — > n + 1 cloning of qubits (see Section Ull fjfl . 



Optimal 1 — > 1 + 1 + 1 cloning of qubits is analysed (Sectio n! H "Fiji . 

A probabilistic implementation of some such asymmetric machines is proposed. (Section IIV|) . 



II. OPTIMAL ASYMMETRIC CLONING MAPS 



Before going into the details of cloning maps, let us sketch the idea behind the subsequent analysis. Let us consider 
the simplest case of 1 — > 2 universal cloning of qubits. Any one-qubit state can be represented by a unit vector n on 
the surface of a sphere, the so-called Bloch vector on the Bloch sphere. The effect of an optimal symmetric cloning 
machine is to produce two output clones with reduced Bloch vectors 7711 (rj < 1). That is, the cloning machine merely 
shrinks the input Bloch vector, but doesn't affect its orientation. This "isotropy" property results from the fact that no 
state is preferred by an optimal universal cloning machine. Therefore the quality of the cloning process is completely 
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characterised by the shrinking factor 77. As demonstrated in pfl, a similar argument applies for N — > M cloning of 
qudits: the quality of the clones is fully characterised by some quantity u, and the whole problem of optimal cloning 
amounts to extremise this quantity. As we will observe below, this argument still applies for asymmetric cloning. 
Consider N — > M cloning and let Ma + Mb — M denote a partition of M. If we require Ma clones to be of the same 
quality, and Mb clones to be of the same quality, the quality of the first (respectively second) set of clones is fully 
characterised by some quantity oj a (respectively uo B ). Therefore, from an algebraic point of view, the main problem 
posed by optimal asymmetric cloning is to find some tight relation between uj a and lo b . 

The object of this section is threefold: (i) to cast the cloning problem into purely algebraic terms, (ii) to give a 
quantitative relation between the qualities of the two sets of clones, (iii) to show how to get optimal cloning machines. 



A. Cloning maps and figures of merit 

An asymmetric cloning device would be an operation taking N identical replicas of an unknown state and output- 
ing M = Ma + Mb approximate clones. In Schrodinger picture, this operation is described by a trace-preserving 
completely positive map fl6| 

T' : T,(M® N ) -> S(J{® M ). 

3i = C d denotes the Hilbert space of a single qudit, IK®^ its iV-fold tensor product, and I B('K® N ) denotes the space 
of bounded operators over 3{.® N . T" makes states evolve and leaves operators invariant . Equivalently, cloning can be 
described in Heisenberg picture by a unital (T(l) = 1) completely positive map T : t B('K® M ) — > 23 (3i® N ). T makes 
operators evolve and states are now left invariant. The two pictures are equivalent and are related by the identity: 

tr( P T(0)) - tr(OT'(p)), (1) 

for all operator G 23(iK® M ) and for all density operator p e r B(^ N ). Note that since all input systems are of the 
form P® N = \i/i)(ip\^ N , IV*) € ^ we can without loss of generality consider the input Hilbert space to be only the 
symmetric subspace of !K® N : 3{® N . So, cloning can be described by a unital cp-map 

T : 3(M <8M ) ■B(0i® Ar ). (2) 

There are essentially two options to quantify the quality of the clones. The first one consists in considering global 
figures of merit. That is if ip denotes the input state to clone, one can consider the quantity 

F A (T) = mbx^ M -\Tr^u B T'(P^ N )\^ M -) 

= min^, (ip® N I (T(P^ Ma ® l^)\^ N ), (3) 

to quantify the first set of Ma clones, where Tr ^®m b denotes the partial trace over the second set of clones. A similar 
expression holds for the global fidelity of the second set of clones, F^UT). 
Alternatively, one can consider single-copy fidelities: 

F A {T) = min^ eM min 1 < fc < M ,(V'|Tr fc r'(|^ Ar )(^ Ar |)^) 

= min^ eM , fe (V^ JV |(T(l^ fc - 1 ® IVX^I ® l^A+M B - k) ^N^ (4) 

for the first set of Ma clones. In this expression, Tr & denotes the partial trace over all clones but the fc-th in the first 
set. A similar expression holds for the single-clone fidelity of the second set of clones, F B (T). 

Even in the symmetric case, it is highly non-trivial to prove that a quantum cloning machine which is optimal for 
one figure of merit is optimal for the other 0. In this work we choose to consider single-clone fidelities because it 
is more relevant than the global fidelities when one considers connections of cloning with other tasks, such as state 
estimation (see Section UlI C|) . or quantum cryptography |l3j . 

Considering the first set of clones, it is obvious that F A (T) lies between the fidelity obtained when preparing the 
clones in a random state, and the fidelity of the optimal N — > Ma symmetric cloning machine F A m (N, Ma)- The 
problem of finding optimal N — > Ma + Mb asymmetric cloning can now be clearly formulated: To maximise F B (T) for 
a given value of F A (T) between 1 /d (the fidelity of a machine producing outputs in a random state) and F A m (N, Ma) 
(the fidelity of an optimal N — > Ma symmetric machine). 

Before proceeding any further, let us fix some notations. These notations are similar to those in the paper of Keyl 
and Werner 0, so as to make the connection between their work and ours as transparent as possible. U(d) will denote 
the (compact Lie) group of unitary d x d matrices, and SU(cZ) will denote the subgroup of U(g0 of elements whose 
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determinant equals 1. The representations of U(d) will be denoted 7r Q . Of particular importance are: (i) the natural 
representation, i.e. the representation given by the elements of the group themselves, which we will simply denote as 
7r, (ii) the TV-fold tensor product of ir, ir® N and (iii) the irreducible representation given by the restriction of tt® n to 
the symmetric subspace J{® N of 3{® N : tt^. Thus, if V>o <= C d denotes some fiducial state, the set of input states to 
clone reads 

{ir+(u)\i) )® N : uSU(d)}. 

In the special case d — 2, as is well known, the irreducible representations are labelled by half-integer positive numbers. 
Accordingly, these irreducible representations will be denoted 7ro, tti/2, etc. u(d) will denote the Lie algebra associated 
to U(d), i.e. u(d) is the Lie algebra of antihcrmitian matrices, and su(d), the algebra consisting of traceless elements 
of u(d). The irreducible representations of u(e?) will be denoted dir a . Let X denote an element of u(d) and let 
{e tx G U(e?) : t e R} denote the associated one-parameter subgroup of \J(d). We have 

n a (e tx )=e td ^ x \ 

Also note that the representation of dn® 1 ^ can be expressed very simply from the natural representation dir: 

M 
i=l 



B. A natural conjecture 

The optimal symmetric cloning map reads (in Schrodinger picture) @: 

T s ' ym : V(Xf N ) -> 23(W® M ) : p® N -> ffl <? M (p® w <g> 1* M - N )S U , (5) 

where d[N] = dimMf N (the constant d[N]/d[M] ensures that the map is trace-preserving). Sm is the projector 
onto the symmetric subspace IK? M . The interpretation of this map is quite intuitive: M — N states containing no 
(quantum) information are appended to the input, and the resulting state is symmetrised. 

Note that the representation 7r® M decomposes as 7r® M = n M 7r rcst , where 7r rost is some representation containing 
no representation equivalent to tt m and acting on some space !H res t E3- Accordingly, we have ]® M = Sm + lrestj 
where l rcs t denotes the identity over 3~C rcst . Thus, what the cp-map © suggests is that optimality is achieved by 
keeping only the component of the decomposition corresponding to the symmetric subspace. 

In the case of asymmetric cloning, we have to consider the decomposition 7t® Ma ® ?t®' Wb w (^tu ® % B ) © 
( n M A ® "rest') © (TTrost ® 7r^ B ) (7T res t ® 7r rcs f)- According to this decomposition, we have 1® M = 1® Ma ® 1® Mb = 

Sm a ® SjVfis + -Sm a ® lrcst' + lrcst <8> 5m b + lrcst ® lrest'- 

Our conjecture is that only the piece t^ Ma <8> 7r} /e should be considered in this decomposition. More precisely, 

let tt Ma <8 n M B * ©t 71 ^ denote the decomposition of t^m a ® "'Mb irreducible components, and let -E 7 denote 
the projector associated to the irreducible component 7r 7 . We conjecture that optimal asymmetric cloning machines 
should be of the form: 

T^ ym : %CKf N ) £(5{® M ) : p®^ -> F*(p® w 1®"-^)^ (6) 

where V is a linear combination of projectors i? 7 . This conjecture is supported by all asymmetric machines we have 
considered. 



C. How to get optimal asymmetric cloning maps 



We now describe a general recipe to get optimal asymmetric cloning machines. This section summarises the results 
obtained in Appendix 1X1 N — > Ma + Mb cloning of qudits is achieved by a cp-map T : S(J{' glM ) — > CB(J{^), which 
decomposes as 

T= E E r( ai ,a 2 ,/3)T( ai ,a 2 ,/3), (7) 

ai£Ll(7r® M A) Q 2 eD(ir® M s) /3GDiv(7rc, 1 «'7r Q2 ) 
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where we define D(ttz) = {a : n a C ttz} and Djy(irz) = {a ■ C nz ® 7r a }. The quantities r(«i, a%, (3) satisfy 

r(ai,a 2 ,f3) > 0, r(ai, a 2 , /3) = 1. (8) 

ai,a 2 ,/3 

The (single-clone) fidelities are essentially fixed by two quantities, to (T) and ui B (T), analogous to the shinrking 
factor r) discussed at the beginning of Sect [HI We have 

F A {T) = 1(1 + ^ A (T)(d-l)), (9) 
F B (T) = l(l + JLu B {T){d-l)). (10) 

The quantities u (T) and uj b (T) are decomposed according to the convex decomposition of T into irreducible 
summands: 

^( T ) E E E r( ai ,a 2 ,/3)^(T( ai ,a a ,/3)), (11) 

" B ( T ) E E E Kai^a.^a^CTCai.aa^)). (12) 

The quantities w" 4 (T(ai, a 2 , /3) are given by 

where C 2 (7r ) denotes the Casimir number associated to the irreducible representation 7r a . The quantities A a satisfy 

Ei a °i 2 = 1 - ( 14 ) 



Similarly, 



B 

UJ a 1 ,a 2 ,P 



1 I |2/1 , ^(^J ~ C 2 (7T b ) 

9 E NU+ 7777+^ }. ( 15 ) 



2 i£D(, 01 ®, ? ) CaC^Jv) 



with £ 6 N 2 = l. 

The quantities A a and /if, are related by 

~r a 

T4 and Wj, are intertwining operators defined as follows. Let X denote an auxiliary space supporting the representation 
[3. Consider the decomposition theory of tt a i ® tt b 2 ® 7r^ w 7r^ © 7r rcst , where 7r rest contains no copy of 7rjy. V a and 
Wb are defined as the unique isometries such that 

«®*a X )V a = V a n+, (17) 

= W b TT+. (18) 

It is now clear that the problem of finding optimal cloning machines is a constrained optimisation problem: We have 
to maximize ui B (T) for a fixed value of F A (T) taking into account the constraint lj5|)- (|14|) - (|16fl . 
In the next section, we apply the recipe just described to treat some concrete examples. 



III. SOME ASYMMETRIC CLONING MACHINES 
A. The simplest example: 1 — > 1 + 1 cloning of qubits 



Although 1 — > 1 + 1 cloning machines of qubits have been extensively studied 01 > it is instructive to revisit them 
in order to illustrate as simply as possible the foregoing analysis. 
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Adopting the standard convention of denoting irreducible representations of SU(2) by half-integer numbers, we have 
here D{tt® Ma ) = D{t:® Mb ) = {1/2}, and £>i(l/2(g) 1/2) = {/3 : 1/2 C 1/2 <gi 1/2 <g) /3} = {3/2,1/2}, where use was 
made of the Clebsch-Gordan series 

for the decomposition of the tensor product of two irreducible representations. Accordingly, 

lo a {T) = r(l/2, 1/2, 3/2)^(T(l/2, 1/2, 3/2)) + r(l/2, 1/2, l/2)^(T(l/2, 1/2, 1/2)), (19) 
u B (T) = r(l/2, 1/2, 3/2)^ s (T(l/2, 1/2, 3/2)) + r(l/2, 1/2, l/2)w B (T(l/2, 1/2, 1/2)). (20) 

With C 2 (j) = j(j + 1), we have: 

^(T(l/2, 1/2, 3/2)) = c B (T(l/2, 1/2, 3/2)) = 1(1 + ^^p) ) = 

So the map T(l/2, 1/2, 3/2) is useless for cloning, since the quantities us A and uj b it yields are worse than for a map 
which would consit of preparing the clones in a random state (w^ ndom = uj b „dom = 0). Let us consider the other 
map, T(l/2, 1/2, 1/2). We have 

^(rd/2,1/2,1/2)) = iflAol'P + Ca(1/ ^ <1/2> ) + fePd + C2(1 ^ (1 ' ) 

= |(2|A | 2 - fl^il^J- ( 21 ) 

Similarly, 

c B (T(l/2, 1/2, 1/2)) = l(2| Al0 | 2 - l^l 3 ). (22) 

Let us now work out the relation between the coefficients Ao,Ai and /zo,/xi. In turn, this relation will give us the 
trade-off between the fidelities of the clone A and the clone B. This relation involves four isometric interwiners: 
Vo , V\ , Wq and W\ . Explicitly, we have 

V = |(l/2-V*)l/2,m}(l/2,m| = Cf™^ 'c[l{ ) 2u){1/2v) \l/2, S ) A \l/2,u) B \l/2,v) x (l/2,m\, (23) 
V X = |(l/2 A l s *)l/2,m)(l/2,m| = C^l^ C^ 1 ' 2 ^ |l/2, S ) A |l/2, u) B |l/2, t,)x(l/2, m|, (24) 
W = |(l/2 S AX )l/2,m)(l/2,m| = Cf^ * ^g^^lA u) A |l/2, S ) B |l/2, u) x <l/2, m|, (25) 

Wi = |(l/2 B l AX )l/2,TO>(l/2,m| - Cg/^ 1 ^ C^ 1 / 2 ^ CV2w) 11/2, w>^| 1/2, S > s 1 1/2, <l/2, , (26) 

where Cv**^**^ denote Clebsch-Gordan coupling coefficients and sum over repeated indices is understood. In this 
expression, |j, m) denote elements of an orthonormal basis for a spin-j representation of SU(2), and |(jij2)j, m) denote 
elements of a basis for the irreducible spin-j representation contained in j% j'2. From Eas l|23l) - (|26|l . one can verify 
that trWfVo = \/3 and trWj*Vi = -1. We can, without loss of optimality, assume that Ao,Ai,/io and £ii are real. 
Using |A | 2 + |Ax| 2 = 1 and |/i | 2 + | Ml | 2 = 1. We get 

^(1/2,1/2,1/2) - 1-^A 2 , (27) 

c B (l/2, 1/2, 1/2) = 1 - ifi = 1 - | (^/T^A 2 - ^Ar) 2 . (28) 

The corresponding fidelities are 

F A = (l+^- 4 (l/2,l/2,l/2))/2, (29) 

F B = (l+w s (l/2,l/2,l/2))/2. (30) 

On Fig^ we have plotted the locus {{F (Ai), F^Ai))}. One readily checks that this locus of couples of fidelities 
correspond with the results in |12| . 
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B. n — > n + 1 cloning of qubits 

We now solve the 2 — ^2+1 case and give the n — ► n + 1 conjectured optimal fidelities. The obtained figures of 
merit agree with the expected values at the limiting points. 

Let us start with the 2^2 + 1 derivation. First, let us observe that D(nf^ 2 ) — {ttq,tti}, D(tt 1 / 2 ) — {^1/2}^ 

£ ) 2(ti"o ® 7r i/2) = { 7r i/2i 7r 3/2}i D2(ni <X> 71-1/2) = {tti/2, 7T3/2, 7i"5/ 2 }- Therefore, according to Sect III Cl the optimal map 
we are looking for can be decomposed as a convex sum of five maps as: 

T = r(0, 1/2, 1/2)T(0, 1/2, 1/2) + r(0, 1/2, 3/2)T(0, 1/2, 3/2) + r(l, 1/2, 1/2)T(1, 1/2, 1/2) + (31) 
r(l, 1/2, 3/2)T(l, 1/2, 3/2) + r(l, 1/2, 5/2)T(l, 1/2, 5/2). (32) 

Consider the map T(0, 1/2, 1/2). This map is characterised by an intertwining operator V satisfying 

VlTl = (ll A <gi 7rf /2 <g> tt^ 2 )F. (33) 
There exists a (Clebsch-Gordan) matrix C such that 7rf /2 <8>7r^ 2 = C* (ir^ x ® 7rf X )C. With V = CV, we thus have 

Vtti = (ir A ® (Trf* 8 Trf x ) )V. (34) 
From Ea. (|13|l . one then sees that 

1/2, 1/2) = i(i + C2(7r ° ) 7? ( " l) ) = °- ( 35 ) 

One can also compute that 

lu b {0, 1/2,1/2) = 1/2. (36) 

Performing the same analysis for the map T(0, 1/2, 3/2), one finds that w A (0, 1/2, 3/2) = and w B (0, 1/2, 3/2) = 
— 1/4. Since a trivial cp-map that merely prepares clones in a random state achieves lu a = oj b = 0, we see that this 
latter map is useless for cloning. The map T(l, 1/2, 5/2) also turns out to be useless because uj a (1, 1/2, 5/2) = —1/2 
and uj b (1, 1/2, 5/2) = —1/4. Consider now the map T(l, 1/2, 1/2). The intertwiner that characterises this map (and 
that we will again denote V) satisfies 

Ftti = (tt a ® 7rf /2 <g> 7rf /2 )K (37) 

Again there exists a (Clebsch-Gordan) matrix C such that V' — CV satisifies 

V'm = {ir A ® {n BX Trf X ))F'. (38) 

The space of solutions for Ea. i|38[) is 2-dimensional: V' is a linear combination of an intertwiner Vq between m and 
tt a (g) tt bx , and an intertwiner V{ between m and 7rj* ® 7rf x . Accordingly, one finds that 

^(M/2,1/2) = |A(1, 1/2, 1/2, 0)| 2 |(1 + Ca(1 ^ a(0) ) + |A(1, 1/2, 1/2, + ) (39) 

= |A(l,l/2,l/2,0)| 2 + i|A(l,l/2,l/2,l)| 2 , (40) 

where |A(1, 1/2, 1/2, 0)| 2 + |A(1, 1/2, 1/2, 1)| 2 = 1 (because V is an isometry). Considering the second set of clones, 
on finds that 

lu b (1, 1/2, 1/2) = 1/2, 1/2, 1/2)| 2 - i| M (l, 1/2, 1/2, 3/2)| 2 , (41) 

where \fi(l, 1/2, 1/2, 1/2)| 2 + 1/2, 1/2, 3/2)| 2 = 1. A similar analysis of the map T(l, 1/2,3/2) shows that 

W A (1, 1/2,3/2) - |(|A(1, 1/2,3/2, 1)| 2 - |A(1, 1/2, 3/2, 2)| 2 ), (42) 
and where |A(1, 1/2, 3/2, 1)| 2 + |A(1, 1/2, 3/2, 2)| 2 = 1, and 

c B (l, 1/2, 3/2) = i|/i(l, 1/2, 3/2, 1/2)| 2 - 1/2, 3/2, 3/2)| 2 , (43) 
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where |/i(l, 1/2, 3/2, 1/2)| 2 + 1/2, 3/2, 3/2) | 2 = 1. Extremising numerically, we have found that optimal asym- 
metric cloning machines are such that r(l, 1/2, 1/2) = 1. Adopting lighter notations, the fidelities can be written as 



F A = 1 



x 

T 



where x 2 + y 2 = 1. Accodring to Sect III Cl optimal machines should be of the form (Schrdinger picture) 



T op t : B{H\ 



®3\ . (gin 



(a* £3/2 + fo*^i/ 2 )(^ 2 ® l)(a£7 3/2 + 6^1/2), 



(44) 



(45) 



where E3/2 and -E1/2 are the projectors onto the irreducible subspaces obtained from the decomposition -k\ ® 7Ti rs 
T3/2 © 7r i/2- Also, 2(2a 2 + & 2 )/3 = 1 because T opt has to be trace-preserving. Note that the fidelities for the limiting 
cases are recovered, that is 



Fx 
F x 



Fo = 1 



F 2 
F 2 
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F3 - -F sym (2 — > 3) — — 
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F sym (l - 2) = - 



^ = 1. 



(46) 



The previous computation strongly supports that the conjectured machines found by combining linearly the pro- 
jectors Er n _x)/2 an d Fr n _x)/2 coming from the decomposition 7r„/ 2 ® 7r i/2 ~ 7r (n-i)/2 © 7r (n+i)/2- The corresponding 
fidelities are given by 



Fx — F2 



n(n + 2) 



F n +x — 1 



1 



n + 2 



(47) 



with x 2 + y 2 = 1. 

These fidelities are depicted in figure ^for n = 1 , . . . , 3, where only the relevant part of the curve is shown. The 
extreme cases are now 



Fx 
Fx 

Fx 



Fo 



F n 



1 



F, 



n+X 



F 2 = ... = F n+ x = F s (n -> n + 1) 



Fo 



F n = F s (n - 1 -> n) 



1 

2 ^ 
_ n 2 + 3n + 1 
~ n 2 + 3n + 2 

n 2 + n — 1 
n(n + 1) 



V = 



! u 2(n+ 1 

= 1 y 



' n + 2 
2(n+ 1) 



(48) 



When n — > 00, an unlimited number of copies of the initial state are available. Then, it is possible to completely 
determine it and prepare a new identical copy. Then F n +x — > 1 as expected. 




0.6 0.7 0.8 

Fidelity for the first N clones 



FIG. 1: Fidelities for the n — > n + 1 cloning machines. The curves for n = 1, 2 are known to be optimal. For n > 2 the curve is 
conjectured to be optimal. Dots on the curves correspond to symmetric machines. 
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C. 1 — > 1 + n cloning of qubits 

The method described in Section UTO has been applied to the case of qubits where one set, A, consists of one clone, 
and the other set, B, consists of n clones (n > 2). 

The optimal fidelities are derived in Appendix iBl and are given by 

F A = l-\y\ F B = l + Ltf + y /^+2jxy), (49) 

3 2 3n 

where x 2 + y 2 = 1 . 

One can check that these fidelities correspond to those obtained using the conjectured form of the optimal cloning 
map. 

Let us discuss these expressions. First of all, imposing F A = 1 implies F B = 1/2. This fact is consistent with the 
idea that in order to prepare a perfect clone, one has to take it from the input and not let it interact with any system 
0. Then, no quantum information is available to prepare the N supplementary clones and the best one can do is to 
prepare the N qubits in a completely random state, thus achieving a fidelity 1/2. 

Second, if one requires the clones B to have the fidelity of an optimal symmetric 1 — > n cloning machine, namely 

F B = ^jt±, one finds that y = y^j^, which gives F A = 3 ^"+^ . Interestingly this fidelity is larger than 1/2: 

in order to produce n optimal clones from one input, not all quantum information need be used and some quantum 
information remains to prepare a "good" (n + l)-th clone. 

Let us now turn to the case of large n. It is well known that there are deep connections between cloning and 
state estimation [fil lift fisj. In particular, for universal symmetric cloning, it appears that there is a correspondence 
between n — > oo cloning machines and state estimation devices [l9| . Such a relation still holds in the asymmetric case. 
Following the lines of [19j, one finds that, in the limit n — > oo, asymmetric 1 — > 1 + n cloning machines interpolate 
between (trivial) machines leaving the quantum system unchanged, and a measuring device estimating destructively 
the input state. In the limit n — > oo, Eqs. H49|) become: 

F A = 1 - \y 2 , F mcas = i + iyv/T^, (50) 

where only the case < y < 1/V2 should be considered. One readily checks that the two extreme cases are found: (i) 
When F A — 1, one finds _F meas = 1/2, which translates the fact that no information can be gained if the input 
state is unperturbed, (ii) The maximum value of F mcaB is 2/3, which is consistent with [2(j. In that case, of course, 
F A = 2/3 too. Between these two cases, the relations (|ST))) express the trade-off between the acquisition of knowledge 
about the state of a quantum system and the disturbance undergone by this system. Actually, such a trade-off had 
been previously studied in 21], in the form of an inequality. So, our machine provides a concrete means to achieve 
measurements saturating this inequality. 



D. 1 — ► 1 + 1 cloning of qudits 

1 — >• 1 + 1 asymmetric machines were first introduced in [12J. The main interest of such machines is that they are 
useful in assessing the security of quantum cryptographic protocols 22]. However, such machines were known to be 
optimal only in the case of qubits. In principle, we could apply the method presented in Section III CI to prove the 
optimality of these. We did not perform such a calculation. Alternatively, one could prove the optimality of these 
asymmetric machines using the isomorphism between CP maps and positive semidefinite operators 23, 24]. Perhaps 
not surprisingly, one finds that optimal cloning machines are of the form |JBJ. Under 7r® 2 , IH® 2 decomposes as a 
-dimensional symmetric subspace JfJ and a -dimensional anti-symmetric subspace !H^~: 

7T® 2 «7T^ ®Tt2, J{^ 2 W^®^, 1® 2 = S 2 + A 2 . (51) 

Let i = . . . d — 1} denote an orthonormal basis of = C d . Clearly, 

1 d-l d-l 

&=2.? (|i>( ' 1 |J>01 + ® |J><Z|) ' M = 2 ? (|i>< *' ® li)01 " |i)01 ® (52) 

i,j—0 

According to Eq.JHJ), in Schrodinger picture, the optimal cloning map is of the form 

T' ovt : S(5C) -» £(J{® 2 ) : p -» (a*S 2 + P*A 2 )(p ® l)(aS 2 + 0A 2 ). (53) 



9 



Since T^ pt should be trace-preserving, we have 

*yV + ^±|/3| a = l. (54) 

Since the map T^ pt is covariant, the fidelity is the same for all input state, and can be calculated with a particular 
state, say |0)(0|. Straightforwardly, we get: 

F A = tr(|0)<0| <g> l)^ pt (|0)<0|) = ^\a\ 2 + + ^(a*/3 + a/?*)- (55) 

F B = tr(l <g> |0)(0|)T o ' pt (|0)<0|) = + ^|/?| 2 - ^i( a */3 + a/?*). (56) 

(57) 

Direct calculations show that these machines correspond to the universal asymmetric cloning machines of qudits 
introduced in |12|. 

E. 1 — > 1 + 1 + 1 cloning of qubits 

We now turn to asymmetric cloning machines with more than two sets of clones. The simplest example of such 
a cloning machine is a 1 -> 1 + 1 + 1 cloning machine of qubits which we shall exhibit now. Following Sect, InTI 
our first task is to determine the representations tt^ satisfying nf? 2 D 1*1/2- O ne directly finds np — 7r , tt\ , 7r 2 . 
Accordingly, optimal 1 — >■ 1 + 1 + 1 cloning maps are of the form 

T = r T + r%Ti + r 2 T 2 , 

where 7> : S(J£f 7 3 2 ) -» B(JC 1/2 ) : A -» V£(j4 ® 1^, with (tt®, 3 <g> Tr^Vg = V^/a, and where r ,n,r 2 > 
0, ro + ri + r 2 = 1. The fidelities for the three clones are given by 

F A = \(l+u A ) F B = l(l + o; fl ) F c = l (1 + wC) 

where 

w A = row + riuj A + r2o; 2 , w s = tqujq + riwf + r 2 w.f , uj c = r wf + nuj + r%uj% . 

Let us work out expressions for u A , u>q , u)q . Let us consider the operator Vq- It satisfies 

(l*l/2 ® ""f/2 ® 7r l/2 ® To^Vb = VbTTi/2. 

Considering the reduction order A ® B ® C ® X — > A (8 -BC (8 X — » ^4 <S> BCX, one finds that Vb is of the form 
Vo = XqEqo + Aii? i, where |A | 2 + |Ai| 2 = 1, and where 

Using Eq. l|13[l . one finds that 
/ , A — 1 4 

W — 1 - -A-l 

Similarly, considering the reduction order A® B ® C ® X ^ B ® AC ® X —> B ® ACX (resp.A ® B ® C ® X ^ 
C <g> BA ® Jf -> C ® BAX). one finds that 

r 4 „ r 4 , 

^o=l- 3M1, W( r = 1 - 3???, 

where |/io| 2 + |mi| 2 = 1, and where |?7o| 2 + |t?i| 2 = L According to Eg. nitty , the following relations hold 

\/3Ao — Ai V3Ao+Ai , . 

Mi = o > = o ■ ( 58 ) 
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With a similar reasoning for the map Ti, one finds that 

<*£ = 1 - |(A? + A?,), wf = 1 - ^(Ai + wp = 1 - | + (59) 

where 

|/* | 2 + lMi| 2 + l£il 2 = i, (60) 

Ao + Ai' _ — Ao + \/2 Ai + Ai' 

Ml = -^- Ml ' = 2 

Ao — Ai' _ — Ao + \/2 Ai — Ai' , . 

7/1 = -7T" m ' = 2 • (61) 

One also easily checks that us^ = = w^F = —1/3, so that the map Ti is useless for optimal cloning, and we can 
choose r2 = 0. In summary, optimal cloning machines are found after maximizing 

A 4A 2 , - 1 4A 2 - 1 

B (A + V3A!) 2 -1 (A + \/2 Ai — Ai<) 2 — 1 

W =ro 3 +ri 3 ' 

c (A - \/3Ai) 2 - 1 (A + V2A 1+ AY) 2 -1 
w = m g +n , (62) 

subject to the normalization, and constraints H58|) . l|61[l . 

Numerical calculations suggest that the optimal solution corresponds to n = 1 (fo = 0). In this case, the optimal 
1 — * 1 + 2 machine given above can be recovered taking Ai' = 0. Interestingly, one finds in that case that some 
quantum information still remains to produce a non-trivial third clone. 

Remark: it is very natural to think of using 1 — > 1 + 1 + 1 cloning machines to perform a simultaneous measurement 
of the three Pauli operators, measuring each Pauli operator at each output of the cloning machine. However such a 
measurement will not be optimal, as has already been demonstrated in |25| using a symmetric 1 — > 3 cloning machine. 



IV. OPTICAL IMPLEMENTATIONS 



We now turn to the issue of implementing some of the machines presented in the previous sections. We will restrict 
ourselves to optical implementations where qubits are represented by polarisation states of photons: N identical qubit 
will be represented by N photons in an identical polarisation mode. In the case of symmetric N — > M cloning, such 
implementations have already been proposed j2(|, and demonstrated experimentally p7|. Let us first briefly recall 
how these N — > M cloning machines work. Let 

|^)-^=(ao^, + /3a^,) w |vac) > (63) 

denote the input state to clone, where |vac) denotes the vacuum state. The labels V, H, s stand for vertical, horizontal 
and signal respectively. Cloning is achieved when a mode prepared in a state (|63|) impinges a crystal where a parametric 
down-conversion (PDC) process can occur. The hamiltonian describing this process is of the form: 

H = 7(4, A* - «L a L) + h. c. , (64) 
where 'i' denotes the idler mode. So the state after the crystal is 

\1> a ) oc e- im {aa^ s + ^J^vac). (65) 

Looking at those cases where there are M photons in the signal mode, one can see that the optimal fidelities for 
the N — > M cloning machine for qubits are obtained. Therefore, the successful realization of the cloning machine is 
conditioned on the number of photons at the output. Note that when M photons are observed in the signal mode, 
N of them came from the initial state and M — N were produced at the crystal, which means that there are as well 
M — N photons in the idler mode. These photons are usually called anti-clones. The total number of photons is 
2M - N. 
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A modification of this scheme was proposed in [28| by Filip, in order to obtain the asymmetric 1 — ^ 1 + 1 machine 
discussed above. His scheme is shown in Fig. |2K- After a successful 1 — » 2 cloning (that is, when all the three detectors 
will click), the two clones are split at the first beam splitter, and one of the clones is combined with the anti-clone at 
a beam-splitter of transmittivity T. Long but straightforward algebra shows that the fidelities at the modes 1 and 2 
give the 1 — > 1 + 1 machine, depending on T. 




FIG. 2: Optical implementation of (a) the 1 — ► 1 + 2 and (b) the 1 — » 1 + 1 + 1 optimal asymmetric cloning machines. 




The natural question is now whether this modification also provides the optimal solution for the more general case 
N — > Ma + Mb- Note that the scheme of |2|| gives all N — > M machines, simply conditioning on different number 
of photons at the input and output signal mode. We denote by Mi the number of photons in mode i. As said above, 
the results of imply that conditioned on the fact that there is one photon at the input signal mode, N = 1, and 
one photon at each output mode, Ma = Mb = M a = 1, the optimal 1 — ► 1 + 1 machine is obtained. It is also easy to 
see that this machine is covariant, so all the calculations can be done taking as initial state (a,y s ) N \v&c). Then, the 
state at the output of the crystal when the total number of photons is 2M — N reads 

j)v,s\j) H ,s\j)v,i\M-N-j) HA . (66) 

The simplest generalization of the 1 — ► 1 + 1 result corresponds to N = 1, Ma — 1 Mb = 2. One can check that 
the evolution of the state l|66|) through the beam-splitters, where N = 1 and M = Ma + Mb = 3, gives the following 
fidelities, depending on T, 

4T*-4T + 7 8T*-4T + 3 

12T2-12T + 9 12T 2 -12T + 9 v ; 

These fidelities are shown in figure |3| Only the relevant part for l/2<T<lis depicted. Note that for T = 1 the 
optimal symmetric machine is recovered, as expected. If the transmittivity decreases, the quality of the first clone 
increases, while the quality of the two clones in mode 2 worsens. When T = 1/2 all the information on mode 2 (and 
3) is lost, and a perfect copy of the initial state is obtained at mode 1. 

How can the missing values be obtained? Note that in the previous expressions, the fidelity for the clone in the 
first mode is always larger than that of the two post-selected photons in mode 2. This suggests a way to find the 
remaining part of the curve of optimal fidelities: one has to reverse the post-selection of photons, that is look at the 
cases where Ma = 2 and Mb = 1. In this way, one expects to reproduce the situation where the fidelity for the two 
clones, now in mode 1, is larger than the fidelity for the single clone, now in mode 2. Repeating the calculations, but 
now for Ma — 2 and Mb = 1, one has 

7T3-4T + 4 3T2-4T + 8 

1 9T 2 -12T+ 12 2 3 9T 2 -12T+12 V ' 

The corresponding curve is also shown in figure |21 In this way the optimal 1 — > 1 + 2 case is completely recovered. 
Indeed, when t = 2/3 it is found that the two photons in mode 1 have fidelities 5/6, while the photon in mode 2 has 
fidelity 5/9, as it should be. Remarkably, the fidelities of Eq. (|67|l are the same as in Eq. I|68f) . changing T into l/T. 
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FIG. 3: Clone fidelities for the optical implementation of the 1 — > 1 + 2 optimal machine. The solid, resp. dashed, line represents 
the case when Ma = 1 and Mb = 2, resp. Ma = 2 and Mb = 1. 



Using the same ideas, we also analyzed the case 2 — > 2 + 1. For Ma = 2 and Mb = 1 one finds, putting N = 2 in 

TO 

(2T-1) 2 (2T-3) 2 
1 2 4(4T 2 -4T + 3) 3 4(4T 2 - 4T + 3) ' V ; 

while for Ma — 1 and Mb = 2 one has 

F = F = 1 {T - 2f F = 1 - ( 3r - 2 ) 2 frnl 

1 2 4(3T 2 -4T + 4) 3 4(3T 2 - 4T + 4) ' K ' 

Note that, again, (f T()|) can be obtained from H69JI if T is replaced by 1/T. These fidelities are depicted in Fig.??, they 
indeed correspond to the optimal solution given above. 

Actually, the case N — > N + 1 can also be computed. The obtained fidelities when Ma — N and Mb = 1 are 

F=F=1 (2T-1) 2 (jVT-TV+1) 2 

1 2 (7V + 2)(27VT 2 -2iVT + 7V + l) 3 (N + 2)(2iVT 2 - 2NT + N + 1) ' 1 ; 

while the expressions for Ma = 1 and Mb = N are again given from these quantities after replacing T by 1/T. One 
can check that the obtained fidelities are identical to Eqs. (|47|l . conjectured to be optimal. 

All the previous results give support to the conjecture that all the TV — > Ma + Mb cloning machines are included in 
Filip's scheme, as it happened for the symmetric case 26] . Unfortunately, this is not the case. Indeed, we've checked 
that this scheme does not provide the optimal solution when N = 2, Ma = 2 and Mb = 2. Therefore, we conjecture 
that this modification of the symmetric cloning machine implementation only works for the cases 1 — > N + 1 and 
N — > N + 1, that is when only two irreducible representations appear in the conjectured optimal solution of 10 • 

This optical scheme can also be adapted to realize the optimal 1 — > 1 + 1 + 1 cloning transformation, see Figl^b). 
Here, the output of a symmetric 1 — > 3 machine is made asymmetric by combining some of the clones with anti- 
clones at two beam splitters, with transmittancc T\ and T2. The obtained fidelities F A > F B > F c , depending 
on T\ and T2, are optimal. The three fidelities are equal when T\ = T2 = 1. Other interesting limiting cases are 
(F A ,F B ,F C ) = (1,1/2,1/2) when T x = 1/2, while taking T 2 = 1 gives Eqs. @3 for the 1 -> 1 + 2 case. This 
construction can easily be generalized further. 



V. CONCLUSIONS 



In summary, we have introduced a new class of quantum cloning machines, which, helps to get a better understanding 
of how quantum information can be distributed unequally between several quantum systems. Mutlipartite asymmetric 
cloning machines have at least two interesting applications. First, some N — > Ma + Mb machines have been proven 
to be a useful tool in assessing the security of some quantum cryptographic protocols 0, Q, • Considering 
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n — > n + 1 machines, we have seen that if one wants to produce n clones from an input with a fidelity which is as high 
as possible, some quantum information still remains to produce a non-trivial n + 1-th clone. Also, we have seen how 
1 — > 1 + n cloning allow, in the limit of large values of n, to study the trade-off between the gain of knowledge about 
the state of a quantum system and the disturbance undergone by this system. We have also demonstrated feasible 
optical implementations of some machines. We have seen that the impossibility of perfect cloning translates in the 
spontaneous emission that unavoidably accompanies stimulated emission. 

Several questions remain open. We here list a few of them. First, it would be very desirable to prove the conjecture 
about the structure of optimal cloning maps (to disprove it would turn even more interesting). Second, it would 
certainly be interesting to find cloning machines optimal with respect to the global fidelities as defined by Eq.lJjjJ, 
instead of single-copy fidelities as we did in this work. Do optimal machines coincide for both figures of merit as in 
the case of symmetric cloning Q? Another interesting problem is to find closed formulas for optimal cloning. We 
have the feeling that this problem will not be solvable with the techniques presented here. The reason is that our 
optimisation requires the computation of a so-called Racah coefficient, for which, to our knowledge, no closed formula 
exists in general. Many interesting questions regarding implementations also remain. It is tempting to believe that all 
N — > Ma + Mb + - ■ ■ cloning machines are only limited by spontaneous emission, and can therefore be implemented by 
splitting clones and anticlones produced by stimulated emission using beam splitters. This question deserves further 
investigation. Finally, it would be interesting to perform an (optical) experiment demonstrating the concepts analysed 
here. 
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APPENDIX A: THE SYMMETRIES OF OPTIMAL CLONING MAPS 



We now use the method of Keyl- Werner and exploit symmetries in order to characterize the sought optimal cp-maps 
T ■ In a sense this subsection can be considered as a summary of their method. There is however a " twist" with 
respect to their analysis, due to the fact that there are now two sets of clones with different fidelities. 

We will work in Heisenberg picture, where states are left unchanged and operators are transformed. For our purpose, 
it is convenient to represent the cp-map we look for with the Stinespring dilation theorem [l6|. This theorem states 
that any cloning map can be written as 

T : %{'K® M ) -» 3(H* N ) : -f V*(0 ® 1 X )V, (Al) 
where 1% denotes the identity over some auxiliary Hilbert space, 3C, and where V : D~C® N — > 1K® M <g> % is an isometry 

(11^11 = IMD- 

The figure of merit @ we have chosen is such that the optimal map can of course be sought amongst covariant 
maps: that is maps T such that VO € r B{'K® M ) and Vu € U(d), 

ir+(u) T(Q) tt+H* = T(ir® M (u) ir® M (u)*). (A2) 

Indeed, since Vit £ U(d), the 'translated' map T u : — > 7r^(u)* T(ir® M (u) 7r® M (u)*) 7r^(u) achieves the same 
fidelity as T, F A (T) = F A {T U ). Thus, 



F A {T) = J duF A (T) = J duF A (T u ) < F A ( J duT u ), 



where du denotes the Haar measure over U{d) \7]{jdu = 1). Similarly, F B (T) < F B (JduT u ). This proves that if 
we find an optimal cp-map T opt , then J duT° pt is covariant and is optimal too. Thus, searching for an optimal map, 
we can restrict ourselves to [/(d)-covariant cp-maps. 

This covariance property is the first symmetry property we will use: it merely states that, for the figure of merit 
we have chosen, no state should be preferred by an optimal cloning machine. Covariance simplifies much the analysis 
because it allows us to use the covariant form of the Stinespring dilation theorem |l6j : 
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Theorem 1 If T : r B('H 1 ) — > 23(Jf 2 ) denotes a cp-map covariant with respect to the representations {U(d) 9 u — > 
7Ti(m) G r S>{'K 1 )} and {{7(d) 9 u — » 7r 2 (u) G 23(3f 2 )}, t/ien T is of the form I^Al)) with the auxiliary space % being the 
carrier space of some representation tx-x ofU(d) and V being an "intertwining" operator: 

(tti(u) ® n x (u))V = Vn 2 (u) Vu g U(d). (A3) 

The second symmetry property to use is of course permutation invariance. Consider constructing a cloning machine 
as follows. We apply the cloning machine described by the cp-map T to our input, and then apply some permutation 
on the output clones of the set A, and some other permutation on the output clones of the set B. Clearly, we expect 
that the performance of the obtained cloning machine will be the same as that of the cloning machine we started 
from, whatever the applied permutations are. Let us formalise this property. Let Sym(Myi) denote the group of 
permutations of M4 objects and let Sym(Myi) 3 p — > 8(p) denote a representation of Sym(Myi) acting on ;H® Ma as: 

6(p) : % miA - % miA : 0! <g> . . . ® 4> Ma -> p( i) ® ■ ® P (Af A )- 
For all admissible cloning map T, the permuted map 

T p : -> TUS^ip) <g> 1® Mb )0(<5(p) <g> 1® Mb )] (A4) 
is a cp-map in its own right and yields clones with the same fidelity as T:F A (T) = F A (T p ) 1 \/p G Sm a - The fact that 

F A (T) = V i^(T„) < F A (— !- V T„) 

pGS Mj4 p£S Mjl 

implies that searching for an optimal map, we can focus on Sym(M^)-invariant cp-maps. Similarly, considering the 
second set of clones, we see that we can focus on Sym(A7s)-invariant cp-maps. 

Consider now the restricted map, T A , obtained upon tracing T over the second set of clones, i.e. 

T A . Ep^M^ ^ <B(ftf N ) : -> T(0 ® 1® Mb ). 

Clearly, T" 4 is a ?7(d)-covariant, Sym(M^)-invariant map, with range in B(5£? ). One can prove that such a map is 
non-degenerate Q, that is there exists a constant w A (T) such that 

T A (dir® MA {X)) = w A (T)dTr+(X), (A5) 

for all X G su(d). Non-degeneracy of T A is a manifestation of the "isotropic" nature of the cloning machine. Indeed 
the fidelity for the clones A reads 

d Ma d 

= \ + ^min^n^+(i^i - \))\r N ) 

= ~(l + ^-u> A (T)(d-l)) (A6) 

The quantity ^w A (T) can be interpreted as the so-called shrinking factor [l^. Clearly, the restricted map T B 

associated with the second set of clones is endowed with the same properties as T A and its cloning quality can also 
be characterised by some shrinking factor Nlo b (T)/Mb- It is because of non-degeneracy that we said in Sect ITT1 that 
the clones of each set are charcterised by a single quantity. 

We now show permutation invariance allows to decompose the sought cp-map as a convex combination of simpler 
cp-maps. Let 



% ®m a = 0q;Kq; x ®m a = ^E a (A7) 



denote the decomposition theory of 7t® Ma . E a is the projector onto "K a . Clearly [E a , tt® Ma (u)] = OV11 e U(<f). 
Hence, by Shur's lemma, T(E a ® 1® m b) \ 8 proportional to the projector onto 'Kj } : T(E a ® 1® M *) = r a S N - So 

T a : S(JC ) -> 2(:k£) : -» r^T^ <g> 1® Mb 0£ q ® l 0Ms ) 
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is a unital cp-map. 

Also, the commutant of (the algebra linearly generated by) all unitaries ir® MA (g) is the algebra linearly generated 
by permutation operators 5(p),p € Sym(Myi). 

Clearly, Vp e Sym(M A ),VO e 'B{'K® Ma ) 1 we have 

T[0(S(p) ® 1® M *)} = T[{6(p) ® 1® Mb )0]. 

Now since each E a is a linear combination of permutation operators and since E\ — E a , we have that T[0(E a <g> 
1® Mb )] = T[(E a <g> 1® Mb )0(£ a <g> 1® Mb )]. Hence the following decomposition holds 

a 

Similarly, we can decompose each T a according to the irreducible representations contained in 7t® Mb . Thus, we get 

T ^ ^ ^ ' T ai ,a 2 T ai ,a 2 , (A8) 

ai£D(Mj) a 2 eD(M B ) 

where Z^M^) = {a\ : ir ai C 7t® Ma }, D(M B ) = {a 2 : 7r a2 C 7t® Mb } and the coefficients r Q1:Q2 are positive constants 
summing to unity. 

It remains to decompose each map T ai , a2 using the covariant form of the Stinespring theorem. The convex 
decomposition of T aitCt2 is the same as the reduction theory of njc into irreducibles. Let 1% — 1/3 denote this 
reduction. V*(l ai <g> l a2 ® 1/3)V = r^l =4> T Q1-Q2 t p : — > r^" V(0 ® 1/3)V* is a unital cp-map and the decomposition 

J , «i,t> 2 = y"V/??ai,a2,/3 (A9) 

holds. In turn, this decomposition induces the decompositions 



X] r «l.«2Xl T '' 3W «l.«2,/3 ( A10 ) 
*1,«2 f3 

Yl r «i,« 2 Xi r /" J «i,a2,i3 ( AU ) 
0:1,0:2 /3 

(A12) 

Relation between the fidelities of the clones. We will now characterise the intertwining operator V and 
see how u A and u B are related. Addressing the first problem requires that we take care of the order in which the 
decomposition of a representation into irreducible components is carried out (this will be clarified below). Addressing 
the second problem requires that we can connect these orders of decomposition with each other. 

Consider a single map T aii0l2i p, and let us solve the equation 

Vn+ (g) = « (g) ® < (g) ® nf (g))V (A13) 

We will consider two manners to reduce ® 7r^ 2 ® 7r?. The first manner first reduces the representation n a2 
(associated with the second set of clones, B), with the representation 7173 (associated to the auxiliary system, X), and 
then the resulting representation ir ai (associated with the first set of clones, A): 

A®B®X ^ A®BX ^ ABX. (A14) 

The second manner is: 

A® B® X -> B®AX -> BAX. (A15) 
Let us consider the first reduction order. Then 

m a 

TT ai ® {7T Q2 «> 7T/3} ~ X ^ T^N.i a © ^rest, 
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where m a denotes the multiplicity of in n ai ® n a , and where 7r rcs t contains no copy of 7r^-. Let us suppose that 
m a < 1 Va. The general case is not more complicated to treat, but this assumption will allow us to adopt lighter 
notations, and at least it holds in all cases exhibited in this paper. Then, we can rewrite the last equivalence as: 

7T ai ® {7T Q2 (g) 7173} W (J) 

TTjV a ® '"'rest , (A16) 

a 

and (up to unitaries) V satisfies 

(0^, o (5)®^rest(3))^ = V7r+(ff). (A17) 

a 

From this relation and Shur's lemma, we infer that there exist coefficients A a such that V — ^2 a X a V a , where V a : 
"H® N — > IH®^ is the unique intertwiner between 7r^ and 7r^- a . One can verify the following properties: 

V a *V b = S ab S N , (A18) 
T ai ,a 2 ,0 is unital =» E |A a | 2 = 1, 

a 

where Sjv denotes the identity over 

Non-degeneracy of T Ql , Q2:( 3 is now expressed as 



dv+{x) = 2 |A a | 2 K*(^ Ql (X) <g> l a )K, X S Bu(d). (A19) 



At this point, it is possible to express ui A as a function of Casimir numbers C2(7r a ), as in pj. We get 

<. Q2l/3 - J + ^{^(^J - E |Aa| 2 C 2 (7r Q )}. (A20) 

In this expression, the sum runs over all a 6 Z?jv(7Tq, 2 ® 7173) = {a : C 7r Q2 ® 7173 ® 7r a }. In the case of qubits 
(SU(2)), where irreducible representations are labelled by positive half-integer numbers 7', we have C^iy) — j(j + 1)- 
Explicit expressions of C2(7r a ) for irreducible representations of SU(<i) can be found in |7lll7|. 

What about the second set of clones? Eq. I|A20|) has been derived following the reduction order l|A14(l , and the fact 
that V = ^2 a X a V a . If instead, we had used the reduction order (|A15|) . we would have found that 

b 

where 7r Ql ® 7173 ~ ® b Tr b , and where W b intertwines 7r^ and the (unique) copy of contained in 7r Q2 ® n b (when any). 
Thus we get 

<, a2 . f 3 = \+ on ) +A C2(n a2 ) - 2 \ti b \ 2 C 2 (n b )}. (A21) 

Z ^ ( ^2(,7Tjv J b 

In this expression, the sum runs over all b € D^r(ir ai ® 7173). All we need now, in order to quantify the trade-off 
between the qualities of the two sets of clones, is a relation between the coefficients A a and the coefficients fi b . It is 
easy to find such a relation: just observe that V = J2 a ^V a = J2 Z ^zW z J2 a KW b V a = J2 Z HzW b *W z = fi b S N - 
Thus 

Ma = E X - teW£V a . (A22) 

+ a 

N.B. The quantity trW£V a is known in representation theory as the Racah coefficient. 
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APPENDIX B: CALCULATIONS RELATED TO 1 — ► 1 + n CLONING OF QUBITS 

We are looking for a map T : 23(J{ ® JC®") -► According to Section GTT3 T decomposes as 

D(jr® n ) is given by the decomposition theory of 7r^,™, which is well-known |l7j : 

^2 = E 

s eJ n i,=i 

In this expression, m s denotes the multiplicity of the representation 7r s , and I n — {1/2, . . . , n/2} when n is odd, and 
In = {0, . . . , n/2} when n is even. Accordingly, 

m s 

T = E E r(s,i s )T(s,i s ), 

where T(s,i s ) : 23(%l/2 ® ^s,0 — > ^Jf^), r(s,i s ) > and X] s i 3 '"(•Ms) = 1- Before decomposing the map T any 
further, we remark that, for fixed s, all maps T(s,i s ) are isomorphic. Therefore, as far as optimality is concerned, 
families of cloning machines with the same values of r s = J^. r(s,i s ) are equivalent. This fact allows to get rid of 
the multiplicities of each s and simply write 

T =E E r(s,P)T(s,l3). (Bl) 

Let us characterise £ ) i(7r 1 / 2 <£> tt s ). We have 7ri/2 C tt±/2 ® tt s 7173 -<=>■ 7173 C 71-®^ X 1 n s- There are three cases 
to consider: Case A: s — 0, which yields -Di(7r 1 / 2 <8> tto) = {0,1}. This case occurs whenever n is even. Case B: 
,s = 1/2; which yields D 1 (ir 1 / 2 ® ^1/2) — {1/2-3/2}. This case occurs whenever n is odd. Case C: s > 1/2, which 
yields £ ) i(7r 1 / 2 <8>7r s ) = {s — l,s, s + 1}. This case occurs whenever n > 1. 

Let us start with case A. So suppose that in the convex decomposition of the cloning map, T, some map T(0) : 
S(JC 1/2 ® Jf ) -> S(JC 1/2 ) appears. T(0) decomposes as T(0) = r(0,0)T(0,0) + r(0, 1)T(0, 1). T(0,0) and T(0, 1) 
have the following structure: 

T(0,0) : 0-> V(0,0)*(O(g)l o M0,0), (B2) 
T(0,1) : O^V(0,l)*(O®li)^(0,l), (B3) 

where (tv* ® Tif <g> 7^)7(0, 0) = V(0, 0)7r 1/2 , and (7if /2 ® tt^ <g> 7tf )V(0, 1) = F(0, l)7r 1/2 . Thus, 

^(1/2,0,0) -|{I + Q(1 ^ W )-I. ^(lA0.0) = l < l + a "g,y ) =0. (B4) 

We see that the map T(0, 1) is useless for cloning. 

The case B is straightforward to treat. Suppose now that in the decomposition of T into irreducible summands, a 
map T(l/2, 1/2) : 'B(Jf 1/2 ® Di 1/2 ) -> T>(M 1/2 ) appears. This map decomposes as T(l/2) = r(l/2, l/2)T(l/2, 1/2) + 
r(l/2,3/2)T(l/2,3/2). The maps T(l/2, 1/2) and T(l/2, 3/2) are exactly those encountered in 1 -> 1 + 1 cloning. 
Thus we see immediately from the results of Section Till Al that: 

cAl/2,1/2) = l-^A(l/2,l/2) 2 , (B6) 

c B (l/2,l/2) = 1 - ^(1/2, 1/2) 2 = 1 - - A(l/2, 1/2) 2 - l\(l/2, 1/2)) 2 , (B7) 

where > |A(l/2, 1/2)| > 1, and that the map T(l/2, 3/2) is useless for cloning. 

We now turn to case C. The convex decomposition of the cloning map T now contains terms T(s) : r B{'Ki/2® !K S ) — ► 
£(^1/12). Each of these maps T(s) decomposes as T(s) = r(s, s - l)T(s, s - 1) + r(s, s)T(s, s) + r(s, s + l)T(s, s + 1). 
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Each map T(s. s — 1) reads 

T(s, s - 1) : - V(s, s - 1)*(0 ® l._i)V(«, a - 1), (B8) 

where 

(7rf /2 ® Trf S tt* JVfa s - 1) = V(«, s - 1)tt 1/2 . (B9) 
There exists a unitary Clebsch-Gordan matrix C such that 

(7rf /2 ® (Trf x © . . .))CV(a, s - 1) = CV(s, s - 1)tt 1/2 . 
We deduce that 

+ (BIO, 

There also exists a unitary Clebsch-Gordan matrix D such that 

(Trf © Trf/f ) DV{s, s - 1) = W(s, s - 1)tt 1/2 , 
if s = 1, whereas Ea. (|B9|) imply that 

(^f ® (^5i/ 2 © Tr s A * /2 )) W(s, s - 1) = W(s, a - 1)tt 1/2 , 
for s > 1. We infer that 

^( S)S -l) = H( s + l). (Bll) 

Let us now consider the maps T(s, s). Each such map reads 

T(s, s) : -» V(s, s)*(G © l a )V(s, s), (B12) 

where 

W% ® *? ® a) = V(s, s)tt 1/2 . (B13) 

Again, there exists some unitary Clebsch-Gordan matrix, which we denote again C, such that 

Wya ® © *f * © ■ ■ 0)CV(*. s ) = s )^i/2- 

From Shur's lemma, V(s,s) decomposes as V(s,s) = A(s, s, 0)V(s, s, 0) + A(s, s, l)V(s, s, 1), where |A(s,s,0)| 2 + 
| A(s, s, 1)| 2 = 1, where V(s, s, 0) intertwines tti/ 2 with the (unique) copy of Tr x / 2 contained in tt^ 2 ® Tr,f x and where 

V(s, s, 1) intertwines 7Ti/ 2 with the (unique) copy of "K\j 2 contained in tt^ 2 ® Trf x . Accordingly, we find that 

■A.,.) - |{i + ac/')-(|A(....o)W(o) + w .,.,i)i-a(i)) ) _ , _ 4 |A(sSi 1)|2 _ (B14) 

A similar reasoning considering the second set of clones gives V(s, s) — fi(s, s,s — 1/2) W(s, s,s — 1/2) + /i(s, s, s + 
l/2)W(s,s,s+l/2), where |/x(s, s, s - 1/2)| 2 + |,u(s, s, s + 1/2) | 2 = 1, where W(s,s,s-l/2) intertwines tt 1/2 with the 
(unique) copy of contained in Trf ® Tr^^y 2 and where Vis, s,s + 1/2) intertwines tti/ 2 with the (unique) copy of 
tti/2 contained in Trf ® tt^^. Accordingly, 

b i v 1 ri , C 2 (s) - Ms, s a - 1/2)| 2 C 2 ( S - 1/2) + jgfo a , g + 1/2)| 2 C 2 ( S + 1/2)) 

w ( s ' s ) = 2 {1 + c^m } ' 

The intertwiners V's and Ws are explicitly given by 

V(s )S) 0) = Cj^f a) C^ v > {sw) \l/2,u) A \ S ,v) B \ S ,w) x (l/2,m\, (B15) 

t/( S , S ,l) = C^f^^ 8 ^ (B16) 

1/2) = C^^C^ (B17) 

W( S , S , S +l/2) = cf 1 7^+ 1/2 ^cf^ (B18) 
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The following relations hold between the A's and the /i' s: 

fi(s,s,s- 1/2) = -{X(s,s,0)tiW{s,s-l/2)*V{s,s,0) + X(s,s,l)tiW(s,s,s-l/2)*V(s,s,l)} (B19) 
H(s,s,s+ 1/2) = -{\(s,s,0)tiW{s,s,s + l/2)*V(s,s,0) + \(s,s,l)tiW(s,s,s + l/2)*V(s,s,l)}. (B20) 
From an explicit calculation(using Mathematica), one gets: 

tvW(s,s,s-l/2)*V(s,s,0) = - trW(s, s, a + l/2)*V(s, s, 1) = 2 J ^ + * ^ , (B21) 



trW(s,s,s-l/2)*V(s,s,l) = tr W(s,s,s + l/2)*V(s, a,0) = W t, 1 - ■ (B22) 



s + 1/2 



From which we find 



w B (s, a) = -(1 - |A(s, s, 0)| 2 ) + ~A(s, s, 0)A(s, s, 1)^ + 1). 

We now turn to the third and last piece: the maps T(s, s + 1) : 23(JCi/2 <E> Jf s ) — > ■B(5£ 1 / 2 )- With a reasoning similar 
to the analysis of the maps T(s, s — 1) and T(s, s), one finds that 

.Ar„ „ , ,a ! ri , C 2 (l/2)-C 2 (l). 



^( SiS+ l) = - {1+ -^ 1/2) - } = -l/3- 



and 



s / ,^ l fi , C 2 ( S )-C 2 ( S + 1/2) 1 2 

So, we see that the maps T(s, s + 1) are useless for cloning. 

Extremisation. Let us first consider the case where n is even. Optimal cloning maps are of the form 

n/2 

T = r(0, 0)T(0, 0) + ^(r(a, s - l)T(s, s - 1) + r(s, s)T(s, s)). 

s=l 

Note that ui A (s, s — 1) = — 1/3 Vs > 1 and that uj b (s, s — 1) strictly increases with s. Thus, for the sake of optimality, 
we can choose r(s, s — 1) = Vs < n/2. Also, one can see that optimal maps can be found for A(s, s, 0), A(s, s, 1) > 
and for r(s,s) = Vs < n/2. Introducing lighter notations; r(0,0) = a, r(n/2,n/2 — 1) = b, r(n/2,n/2) = c, 
A(n/2, n/2,0) = x, and A(n/2,n/2, 1) = y, and extremising, we get Eas. iEE5|) . 
A similar argument holds when n is odd (n > 3) and also leads to Eas. l|49ll . 
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